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ABSTRACT 

Extreme helium stars are very rare low-mass supergiants in a late stage of evolution. They are 
probably contracting to become white dwarfs following a violent phase of evolution which 
caused them to become hydrogen-deficient giants, possibly RCrB stars. Using the latest 
generation of models for spherically expanding stellar atmospheres, we set out to measure 
mass-loss rates for a representative fraction of these stars. We have used high-resolution 
ultraviolet and optical spectra, and ultraviolet, optical and near-infrared photometry from 
a variety of archives. Overall atmospheric parameters have mostly been taken from previ- 
ous analyses and checked for consistency. Mass-loss rates were measured by fitting the P- 
Cygni and asymmetric profiles of C, N and Si ultraviolet resonance lines and lie in the range 
10~ ln — 10 Moyr -1 . These rates follow a Castor-type (M oc L 1 * 5 ) relation marking a 
lower limit for the mass loss from hot stars of all kinds. The mass-loss rates of the studied 
stars also show a strong correlation with their proximity to the Eddington limit. There is no 
firm evidence for variability in the stellar wind, although photospheric pulsations have been 
reported in many cases. 
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1 INTRODUCTION 

Amongst the more unusual groups of stars in the Galaxy is a 
collection of low-mass but nevertheless luminous stars from whose 
surfaces hydrogen is almost completely absent. Instead they are 
composed principally of helium and some 3% carbon (by mass). 
The collection includes RCrB stars, with rough spectral classes F 
and G, extreme helium stars, with cl asses A and B, and some he- 
lium subdwarf O stars djeffervll2008l) . It had long been questioned 
whether these stars are hydrogen-deficient because their outer lay- 
ers have been stripped by a stellar wind, as in the massive Wolf- 
Rayet stars, or diluted by deep mixing of the stellar envelope, as 
in the low-mass Wolf-Rayet c entral stars of plan etary nebulae, or 
by some more exotic process l lSchonbernerl 1 198(3) . Increasing ev- 
idence has accrued that these stars are the result of a merger be- 
tween a helium- and a carbon/oxygen white dwarf. In such stars 
the helium white dwarf is entirely disrupted, wrapping itself around 
its more massive companion, whereupon a helium shell is ignited 
and the star becomes a helium supergiant for a few tens of thou- 
sands of years, before contracting to become a massive C/O white 
dwarf JWebbinklll984liben & Tutukovl l l984l : ISaio & Jeffery|2002l ; 
IPandev et aiTl2006l : IClavton et alj|2007t) . While the arguments for 
the white dwarf merger model are compelling, sufficient questions 
remain that it is important to extract as much physical information 
concerning these stars as possible. 
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It has been known sin ce the early days o f the International 
Ultraviolet Explorer (IUE) dKondo et all 1 19891) that some of the 
hotter (B- and O-type) hydrogen-deficient stars showed significant 
P Cygni-type line profiles in strong ultraviolet resonance lines, indi- 
cating the pres ence of a stellar wind and hence mass loss from the 
stella r surface I Darius et al J 1 19791; [Giddingsl 1 1 98 1| : lHamann et al.l 
1 19821 ; IJefferv et alj Il98fj ; iDudlev & Jeffervl Il992l) . Estimates of 
mass-loss rates of between 10~ n and lO _7 M0yr~ 1 from these 
early studies showed that the current stellar wind is unlikely to be 
responsible for stripping the entire hydrogen envelope. 

However, the study of the winds in these stars is important for 
several other reasons. First, it is well known that a relation exists be- 
tween lumi nosity and mass-loss rate (or, more precisely, wind mo- 
mentum) (Kudr itzki et al J 19991) . so that the latter becomes a proxy 
for the former. This is extremely useful when distances are not 
known. Moreover, models for pre-white dw arf contraction give a 
direct relation between luminosity and mass JSaio & Jeffery 1988, 
l2002h . 

The drawback is that the modelling of a stellar wind requires 
a difficult computation in order to satisfy conditions of radiative 
and statistical equilibrium, and momentum and energy conserva- 
tion. Specifically, the conditions of local thermodynamic and hy- 
drostatic equilibrium no longer hold, and it is necessary to compute 
a self-consistent multi-level multi-atomic model in an expanding 
medium. Fortunately, in the quarter century since this topic was last 
addressed, our ability to perform such calculations has improved 
by orders of magnitude, so that a re-examination of the question 
is timely. Additional data from IUE, the Hubble Space Telescope 
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(HST), the Far Ultraviolet Space Explorer (FUSE), the Tubingen 
Echelle Spectrograph aboard the Orbiting and Retrievable Far and 
Extreme Ultraviolet Spectrometer (ORFEUS/TUES) and several 
ground-based telescopes provide an additional opportunity to im- 
prove on earlier results. In the cases of the hottest (O-type) helium 
stars, the LTE assumption is of little value, so that a study of this 
type provides an opportunity to gain new knowledge about the sur- 
faces of these stars. 

This paper addresses the following questions. Do early-type 
low-mass helium stars show conventional radiatively driven stellar 
winds and, if so, what are the associated mass-loss rates? Can such 
a mass-loss rate be used as a proxy for the luminosity, and hence 
mass, of the star in question? Are the results consistent with the 
predictions of stellar evolution theory? 

Section 2 describes the observational data available for such 
a study and explains which data have been used and why. It 
gives more information about the particular stars examined - 
namely BD+37°442, BD+37°1977, HD 160641, BD-9°4395, 
BD+10°2179, and HD 144941. Section 3 describes the theoretical 
framework for the wind analysis, including the self-consistent wind 
models, the formal solutions and the methods used to fit the obser- 
vations. Section 4 presents the measurements for individual stars, 
and Section 5 concludes with an assessment and interpretation of 
the measurements. 



2 OBSERVATIONAL DATA 

To determine stellar wind properties it is necessary to observe 
the profiles of various spectral lines. Most of the lines of interest 
lie in the ultraviolet and far ultraviolet; therefore high-resolution 
spectroscopy from space is required. To support these data, the sur- 
face temperature and chemical composition of the star also need to 
be known. The former is best described by fitting the overall flux 
distribution, the latter by fine analysis of the high-resolution photo- 
spheric spectrum. For three of the stars studied here, measurements 
of both have already been made; these show some discrepancies 
from one publication to another, which can be attributed to im- 
provements in either the observations or the model atmospheres. 
The following subsections describe the stars selected for study, the 
ultraviolet data available for wind analysis, and other supporting 
ground-based spectroscopy and photometry. 

2.1 Sample Selection 

The primary criterion for selection is that targets should be bright 
enough for high-resolution ultraviolet spectroscopy, and so must 
be hot with low interstellar extinction. Th is limited the origi- 
nal study of winds in extreme helium stars dHamannetaUT l982) 
to HD160641, BD-9°4395, and BD+10°2179. Other extreme 
helium stars with effective temperaturf] T e g > 15 000 K in- 
clude BD+13°3224, LSS3184, LSIV+6°2, LSS5121, DYCen, 
LSS 435 7, LSE78, LSS 9 9, LS IV- 14 109, HD124448 and 
MV Sgr Jjefferv et alj|200lh . While all have been observed with 
IUE, none show evidence of strong PCygni line profiles at low 
resolution. The hottest stars (LSIV+6°2 and LSS 5121) are heav- 
ily reddened and have not yet been observed at high resolution. 



1 We distinguish the effective temperature T c ff obtained from analyses of 
photospheric spectra from the surface temperature T„ conventionally used 
in stellar wind models (§ 3). 



Thus the secondary selection criterion is that suitable data are al- 
ready available. Figure Q] shows a part of the spectrum of those 
stars which satisfy these criteria in the vicinity of the Civ reso- 
nance double at 1548-50A. Strong PCygni profiles are evident in 
the hottest stars (top), weakening towards lower effective temper- 
ature (bottom). For BD+10°2179 and HD144941, the Civ lines 
are asymmetric. For BD+13°3224, they are slightly broadened and 
possibly purely photospheric. For HD168476, only the interstellar 
component is readily identified. 

In addition to these classical extreme helium stars we included 
three closely-relate d stars. HD144941 has a hydrogen abundance 
around 3% by mass faarrison & Jeffervl 19971) , and is carbon weak. 
It has been linked to low-luminosity extreme helium stars like 
BD+13°3224, although the nature of the connectio n is unclear. 
The U V spectrum shows evidence of a stellar wind (Jeffery et al. 
fl986h . 

It is not known whether the two helium sdO stars BD+37°442 
and BD+37° 1977 represent very hot extreme helium stars, al- 
though they are known to be hydrogen-deficient ( Wol ff et al.l 19741) 
and very luminous relative to other helium-rich sdO stars. No com- 
plete quantitative analysis has been formally published for either 
star. A non-LTE stud y of the winds in bot h stars gives mass- loss 
rates jGiddingsIl 198 lh and was reported bv lDarius et al.1 Jl979h . A 
photospheric analysis of BD+37°442 by Husfeld was reported by 
Bauer & Husfel 3 ll995h and shows H and C abundances similar to 
other extreme helium stars. 

Finally, we omitted BD+13°3224, HD124448 and 
HD 168476 from our analysis on the grounds that these data 
either showed little evidence for a stellar wind, or were too noisy. 

2.2 Ultraviolet Spectroscopy 

High-resolution ultraviolet spectra of extreme helium stars have 
been obtained with four spacecraft, namely the International Ul- 
traviolet Explorer (IUE), the Hubble Space Telescope (HST), The 
Far Ultraviolet Explorer (FUSE) and the Orbiting and Retrievable 
Far and Extreme Ultraviolet Spectrometer (ORFEUS II). 

A summary of the data available for extreme helium stars with 
T c ff > 12 000K is given in Table Q] In four cases, spectra were 
obtained in overlapping wavelength intervals by more than one in- 
strument, and found to agree well. Where multiple exposures of a 
given object were obtained with the same instrument, the reduced 
spectra were combined to increase the overall signal-to-noise ratio, 
using the individual exposure times as weights. 

The greatest body of ultraviolet spectroscopy for extreme he- 
lium stars continues to be SWP HIRES images obtained with 
IUE. Despite the larger collecting area of the Hubble Space Tele- 
scope, shorter exposure times mean that STIS images which in- 
clude the key wind lines (C III, C IV, Si IV, N V) have lower signal- 
to-noise than the corresponding merged IUE data. ORFEUS data 
for BD+37° 1977 are of considerably higher S/N than the merged 
IUE spectrum, so these data have been used at wavelengths short- 
ward of 1400A. For BD+37°442, the IUE and ORFEUS data are 
of comparable quality. 

2.3 Optical Spectroscopy 

Blue-optical spectra have been used to ensure that the wind mod- 
els are consistent with the global emergent spectrum and, in par- 
ticular, with the profiles of strong lines formed in the photosphere 
and unaffected by the wind. These data come from a variety of 
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Figure 1. The Civ ultraviolet resonance doublet in all extreme helium stars observed with IUE. The rest position of the doublet is marked. The stars are 
arranged in decreasing effective temperature from top to bottom. All spectra are from SWP HIRES images and are the result of merging the available spectra 
for each star (Table[T}. They have been wavelength-corrected to bring photospheric lines close to laboratory wavelengths. The strong feature around 1526A is 
due to interstellar absorption, the sharp features around 1540, 1555 and 1569A are due to fiducial marks on the camera. 



sources (Table QJ and are also used with caution. Three of the tar- 
gets are known or suspected to be variable; the adopted optical 
spectra of HD160641 and BD— 9° 4395 are direct averages of time- 
series spectra obtained over four and three nights respectively. Both 
stars are pulsating and and show substantial line-profile variabilit y 
on a timescale of hours jjefferv & Heberll992l:l Wright et alJ2 006). 
Variability with P ~ 590s in BD+37°442 jBartolini et aUll982t) 
has not been confirmed. 



2.4 Ultraviolet, Optical and Infrared Spectrophotometry 

Finally, to ensure conservation of total flux and to check the 
interstellar extinction, spectrophotometry from the ultraviolet to 
infrared was collated . UV spectrophotometry was the same as 
used by Ijefferv et alj foOOlh . with the addition of IUE images 



SWP05642, LWR04892 and LWR02805 for BD+37°442. Opti- 
cal photometry (Johnson B and V) were taken from SIMBAD, and 
infrared fluxes were retrieved from 2MASS. 



3 THEORETICAL FRAMEWORK 

For analyzing the stars, we compare the observations with syn- 
thetic spectra calculated with the Potsdam Wolf-Ray et (PoWR) 
model atmospheres code ( Hama nn & Grafenerj 12004 and refer- 
ences therein). This code solves the non-LTE radiative transfer in 
a spherically symmetric expanding atmosphere. Mass-loss rate and 
wind velocity are free parameters of the models. 

Wind inhomogeneities ("clumping") are accounted for in a 
first-order approximation, assuming that optically thin clumps fill 



4 C.S. Jeffery & W.-R. Hamann 

Table 1. High-resolution spectroscopic data used in the current analysis 



Star 


Instrument (Observer) 


Image Number 


Date 




A [A] 


Res. [A] 


BD+37°1977 


IUE SWP HL 


06766 


1979-10-05 


3600 


1 150 


- 1980 






IUE SWP HL 


07248 


1979-11-28 


2643 


1150 


- 1980 






MacDonald 2.7m Coude 








3700 


- 7000 


0.15 


Bu-\-o7 442 


TT TT"; OH7T1 ITT 

IUh SWr HL 


06768 


19/9-10-05 


3600 


1150 


- 1980 








S90Q ^ 


1 QQ^ 1 1 in 


2341 


904 - 


- 1410 






Calar Alto 2.2m FOCES 








Jyau 


— oouu 


0.1 


HFl 1 ^0^41 


UCT QTTS n 9V0 9 F14.0M 


UDOVlfiUlU 


9000 10 91 

ZUUU 1 U iJ 


1 50 


1 145 


- 1729 






ESO 3.6m CASPEC (1) 


mean 


1984-04-03, 1985-04-09 




3900 


-5160 


0.15 




AAT 3.9m UCLES (2) 


mean 


2005-08 






— J lOU 


0.15 


du — y 'ioyo 


TT TF SWP T-TT 
1UD j Wr HL 


13826 


1 1 04 98 


23880 


1150 


-1980 






IUE SWP HL 


30797 


1987-04-17 


24240 


1150 


- 1980 






IUE SWP HL 


30808 


1987-04-19 


23520 


1150 


-1980 






IUE SWP HL 


30814 


1987-04-20 


20460 


1150 


-1980 






IUE SWP HL 


30818 


1987-04-21 


21540 


1150 


- 1980 






ESO 3.6m CASPEC (3) 


mean 


1987-04 




3994 


-4988 


0.15 




ESO 3.6m CASPEC (3) 


mean 


1987-04 




5778 


-6790 


0.15 


BD+10°2179 


HST STIS 0.2X0.2 E140M 


O66V04010 


2001-02-11 


144 


1145 


- 1729 






ESO 3.6m CASPEC (1) 




1985-04-08 




3790 


-4810 


0.15 


HD 144941 


IUE SWP HL 


23961 


1984-09-14 


12600 


1150 


- 1980 






IUE SWP HL 


23962 


1984-09-14 


8940 


1150 


- 1980 






AAT 3.9m UCLES (4) 




1995-04-15 




3900 


-5160 


0.15 



Observer: (1) Heber; (2) Ahmad, Sahi'n & Jeffery; (3) Heber & Jeffery; (4) Harrison & Jeffery 



a volume fraction /v while the interclump space is void. Thus the 
matter density in the clumps is higher by a factor D = /y , com- 
pared to an un-clumped model with the same parameters. D =4 is 
arbitrarily assumed throughout this paper. Nothing is known about 
the wind homogeneities in our program stars. Even for the winds 
from massive stars, the clumping contrast is highly debated (see 
Hamann et al. 2008). Fortunately, this issue has no impact on the 
present study, since we derive the mass-loss rates from resonance 
lines. Since their (mainly line-scattering) opacity scales linearly 
with density, inhomogeneity does not change the average opacity 
as lon g as the individual clumps are optically thin dOskinova et al.l 
120071) . Moreover, since the resonance lines are not strongly sat- 
urated in the observed spectra our program stars, the clumps are 
certainly not optically thick. 

Detailed ions of H, He, C II - C IV, N II - N V and Si III - Si IV 
are taken into account in our models. About 350 non-LTE levels 
and 4500 line transitions are explicitely treated. Because of the un- 
settled questions about how to include dielectronic recombination 
properly, this process is neglected. 

The models account for line blanketing by iron and other iron- 
group elements. About 10 5 energy levels and 10 7 line transitions 
between those levels are taken into account in the approximation of 
the "supetievel" approach. 

As already mentioned, the velocity field must be pre-specified. 
In the subsonic region, v(r) is defined such that a hydrostatic den- 
sity stratification is approached. For the supersonic part we adopt 
the usual /3-law, basically of the form v(r) = t>oo(l — 1/V)' 3 > w i m 
the terminal velocity Voo being a free parameter. The exponent f3 
is chosen between 0.7 and 1.0 such that the wind profiles are re- 
produced best. For some program stars we employ a "double beta 
law", which means that a second term of the same /3-law form is 
added to the velocity, but with a large /3 = 4. Typically, 40% of 
the terminal velocity is attributed to this second /3 term. Its effect is 
that the velocity v(r) does not approach the terminal velocity too 
soon. The "double beta law" makes the blue edges of P-Cygni pro- 
files less steep. This often fits the observation better, and has also 



fou nd theoretical support from hydrodynamically consistent mod- 
els jGriifener & Hamannll2005l) . 

Further input parameters to the PoWR models are the lumi- 
nosity L and the "stellar temperature" T*. The latter refers, via the 
Stefan-Boltzmann law, to the radius i?» which is the inner bound- 
ary of the model atmosphere, by definition located at a Rosseland 
optical depth of 20. In contrast, the usual "effective temperature" 
Teff refers to the radius where the Rosseland optical depth is 2/3. In 
thick expanding atmospheres, e.g. from Wolf-Rayet stars, this can 
be considerably larger than i?» . The winds of our program stars are 
rather weak, and hence the difference between T„ and T e a can be 
neglected here. 

The Doppler velocity vo reflects random motions on small 
scales ("microturbulence"). In stellar winds vo is usually found to 
be of the order of 10% of the terminal wind velocity. Unfortunately 
the PoWR code does not allow to vary vo with the depth in the 
atmosphere. Therefore we apply a smaller vo value in the "formal 
integral" when calculating the photospheric absorption spectrum in 
the optical, and a larger value for the UV spectrum with the wind 
lines. The millions of iron-group lines are always calculated with 
vo = 50 km s" 1 for numerical reasons, which is certainly too high 
for the photosphere and therefore may lead to a slight overestimate 
of the iron-line forest in the UV. 

In addition to Doppler broadening from microturbulence, ra- 
diation damping and pressure broadening of spectral lines are 
taken into account in the formal integral. This is a new feature 
of the PoWR code not yet documented in previous p ublications. 
The br oad ening of H I and He II lines uses tables from lVidal et al.l 
dl973l) and lSchoning & Butlej Jl989cllbll3) respectively. Voigt func- 
tions are used for lin es of neutral helium , with calculation of the 
coefficients following Griem et al. Excited lines from other 

elements are broade ned by the qu adratic Stark effect following the 
approximation from lCowlevI dl97 lh . 

Pressure broadening is missing in the code for lines that are 
not covered by the mentioned sources. We have no pressure broad- 
ening included yet for resonance lines, which is a noticeable de- 
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Figure 2. Spectral energy distribution (SED) for BD+37°1977 from the UV to the near IR. The observations comprise low-resolution IUE data (thin blue 
line), visual photometry (blue blocks with magnitudes) and near-IR photometry from 2MASS. The model continuum (see Table[3]for the parameters) is shown 
by the straight red line, in parts of the spectrum augmented by the synthetic line spectrum. Since the model luminosity was adopted from previous work, the 
SED fit is achieved by adjusting the color excess Eg_v an d the distance modulus DM. 



Table 2. Previously published physical parameters for the program stars. Abundances X are mass fraction in per cent. 
Star 

T e g E B _v log 9 logL v x logM Ref. X H X c X N X Si X Fc Ref. 

[kK] [mag] [cms" 2 ] [L ] [kms" 1 ] [Movr" 1 ] [%] [%] [%] [%] [%] 

BD+37°1977 

41.0 0.00 1 

55.0 0.00 4.0 4.4 2200 -8.2 ... - 8.0 2 

-9.8 3 

BD+37°442 

60.0 0.0 4.0 4.5 4 2.5 0.3 0.08 5 

55.0 0.00 4.0 4.4 2200 -8.2... - 8.0 2 

-9.4 3 

HD 160641 

31.6 0.45 1 0.3 0.5 0.06 6 
34 2.8 7 

4.5 550 -8.2... - 7.2 8 

-10.8 3 

BD-9°4395 

22.7 0.24 2.55 4.1+ 1,9 0.04 1.3 0.09 0.14 0.017 9 

600 -8.5... - 7.7 8 

-10.5 3 

BD+10°2179 

19.5 0.00 2.6 3.7+ 1,10 3 0.13 0.04 0.01 10 

16.9 0.00 2.6 3.5+ 11 2.2 0.08 0.013 0.006 11 

400 -11.0... - 8.9 8 

-14.4 3 

HD144941 

23.2 0.25 3.9 2.7+ 12 1.4 0.005 0.003 0.0025 0.002 12,13 

27.8 0.25 1 

350 -9.7... - 9.0 14 

Notes: f: estimated from T c g , log g using the Af corc — L s hcll relation for He-shell burning stars (Jefferv 1988). 
References: 1: Jefferv et al. (2001) 2: Darius et al. (1979) : Giddings (1981) 3. Dudley & Jefferv (1992) 

4: Husfeld (1987) 5: Bauer & Husfeld (1995) 6: Aller (1954) 

7: Rauch (1996) 8: Hamann et al. (1982) 9: Jefferv & Heber (1992) 

10: Heber (1983) 11: Pandevetal. (2006) 12: Harrison & Jefferv (1997) 

13: Jefferv & Harrison (1997) 14: Jefferv et al. (1986) 
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ficiency in our spectral fits. This somewhat incomplete treatment 
of line broadening is due to the historical fact that the PoWR code 
was developed for stellar wind spectra where pressure broadening 
is negligible. 



4 SPECTRAL ANALYSES 

Since the winds from our program stars are weak, their spectra 
mainly form in the nearly hydrostatic photosphere. Except for the 
strongest lines which are affected by the stellar wind, static, plane- 
parallel non-LTE models therefore provide a reasonable approxi- 
mation. All our program stars have been previously studied with 
static models in more or less detail. We adopt the parameters from 
these previous studies, at least as a first guess. This refers to the 
effective temperature T e jj, the surface gravity log g, the stellar lu- 
minosity L, and the chemical composition (see Tabled- We check 
the fit of the photospheric absorption line spectra with our models, 
and keep the parameters from previous work whenever possible, 
but also encounter the need to adjust parameters in a few cases (see 
below). 



4.1 BD+37°1977 

Our fit procedure with the PoWR models is decribed now by the 
example of the first program star, BD+37° 1977. The observed and 
computed spectral energy distribution (SED) is shown in Fig.[2] 
Since the model luminosity was adopted from previous work, the 
SED fit is achieved by adjusting the color excess Eb-v and the 
distance modulus DM to the values given in Table|3] 

The model SED depends of course also on the stellar pa- 
rameters. Hence we have to look at the line spectrum as well, 
and the whole procedure is in fact iterated. The lower panel of 
Fig.|3]shows the optical spectrum, which is dominated by the lines 
from neutral and ionized helium. Getting the balance between the 
He I and the He II lines right constrains the effective temperature. 
Previous (static) analyses give slightly discrepant values for T e g 
(41 kK and 55 kK, respectively, see Table[2j. We find a good fit with 
T e g = 48 kK. The pressure-broadened line wings are reasonably 
well reproduced with log g = 4.0 from the previous work. 

Unfortunately, BD+37°1977 was not analyzed yet for the 
chemical abundances. However, BD+37°442 can be considered as 
a spectroscopic twin, and the composition of the latt er has been 
studied with static models bv lBauer & Husfekj dl995h . Therefore 
we adopt the abundances from that paper. The iron mass fraction 
we set to half solar. With this composition, the simulated "iron for- 
est" in the UV still appears a bit stronger than observed, as does 
the Si IV resonance doublet. But this problem can be attributed to 
the effect of microturbulence. The UV range has been calculated 
with a microturbulent velocity of 100km s^ 1 , which is adequate 
for the P Cygni profiles from N V and C IV that are formed in the 
very fast wind of this star. For the photospheric lines this value is 
far too high, but the code allows only for a uniform value. The opti- 
cal range is calculated with a microturbulence of 30kms _1 which 
may be adequate for the photosphere. 

Finally we turn to the wind lines in the spectrum of 
BD+37° 1977, which are the P Cygni profiles of the N V and the 
C IV resonance doublets shown in the upper panel of Fig.[3] or in 
more detail in Fig.|4] The terminal wind velocity Vrx is adjusted to 
fit the width of the P Cygni profiles, while the mass-loss rate M 
influences their overall strength. The parameters for the best fit are 
included in Table[3] The two wind-formed lines available are close 



to saturation. Hence the error margin of the mass-loss rate (usually 
about a factor of two) is asymmetric in this case, i.e. much higher 
rates cannot be safely excluded. One should also keep in mind that 
any error in the nitrogen and carbon abundances directly propagates 
into our empirical mass-loss rate. 

Note that in this semi-empirical approach, the mass-loss rate 
and wind velocity are free parameters. Hydrodynamicalfy consis- 
tent solutions that agree with observations have been constructe d 
for, e.g. massive O stars (for a r ecent review see Puis et al. 2008), 
and some Wolf-Rayet subtypes jGriifener & Ham ann 2005, 120081) 
under the paradigm of radiation-driven winds. Problems with such 
models are encountered not only for very strong m ass-loss, but also 
for ve ry thin winds (the "weak wind problem", cf. lMarcolino et all 
120091) . A first indication if the empirically found mass-loss is plau- 
sibly driven by radiation pressure is given by the work ratio, de- 
fined as the total momentum intercepted by the wind from the ra- 
diation field over the mechanical momentum carried by the stellar 
wind after escaping from the star {M v^). A computed work ratio 
somewhat below unity, say about 0.5, is expected because of the 
incompleteness of opacities in the model calculations. Note how- 
ever that for being fully hydrodynamically consistent the momen- 
tum balance must be fulfilled not only globally, but at each radial 
point. In the case of BD+37° 1977, our final model has a work ratio 
of 0.60, which means that the parameters are very plausible for a 
radiation-driven wind. 



4.2 BD+37°442 

BD+37°442 is a spectroscopic twin of BD+37°1977, although 
previous analyses assigned a slightly higher effective temperature 
to the former. Within the uncertainties, we can adopt the same best- 
fit model for both objects. BD+37°442 shows more interstellar 
redding, but still appears brighter than BD+37° 1977 which must 
be attributed to a smaller distance if the intrinsic luminosities are 
the same (see TableO. As mentioned above, we c an take the chem- 
ical co mposition from a photospheric analysis by Bauer & Husfeld 
jl995l) . The complete SED and line-fit plots are provided as on- 
line material, while only the zoomed-in wind lines are shown in the 
main paper (Fig.|5j. The different observations for the Nv profile 
(blue: ORFEUS, green: IUE) demonstrate the problems of back- 
ground subtraction in high-resolution IUE data. 



4.3 HD160641 (V2076 Oph) 

The helium absorption lines in the optical spectrum are well repro- 
duced with an effective temperature of 35.5 kK, slightly higher than 
found in previous work. The o nly abundan ce analysis published to 
date is a "quick look-see' ' by lAUerl Jl954h . For the present study 
we set the carbon abundance to 1 % (mass fraction), which leads to 
a reasonable reproduction of the photospheric C III and C IV lines. 
Such a carbon abundance of a few times solar seems to be typi- 
cal for extreme helium stars. For the nitrogen abundance the solar 
value also produces too weak photospheric lines, while five times 
solar (0.3% mass fraction) gives a reasonable fit. The Si IV line at 
4089 A forces us to adopt twice the solar silicon abundance. Solar 
iron perfectly reproduces the UV iron forest, and hydrogen is ap- 
parently absent. The iron forest has also significant impact on the 
SED, which is perfectly fitted by the model. The fit of the SED 
and the UV and blue-optical line spectrum are shown in the on- 
line material (Fig.[2]l. In the formal integral, different microturbu- 
lent velocities are again used for the UV range with the wind lines 
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Figure 3. Normalized line spectrum for BD+37°1977 in the UV (top) and optical (bottom) range. The optical observation (blue line) has been normalized to 
the continuum "by eye", while the absolutely calibrated UV spectrum was divided by the model continuum. The model spectrum (red dotted) is for the model 
parameters given in Table[3] 



Table 3. Current physical parameters for each program star. Abundances X are mass fraction in per cent. 



Star 


T e ff 


E B -v 


log 9 


logL 


DM 




logM 


X H 


X C 




X S i 


^Fe 




[kK] 


[mag] 


[cm s — 2 ] 


[L©] 


[mag] 


[kms- 1 ] 


[Moyr- 1 ] 


[%] 


[%] 


[%] 


[%] 


[%] 


BD+37°1977 


48.0 


0.03 


4.0+ 


4.4+ 


12.2 


2000 


-8.2 





2.5 


0.3 


0.08 


0.06 


BD+37°442 


48.0 


0.09 


4.0t 


4.4+ 


11.7 


2000 


-8.5 


0+ 


2.5+ 


0.3+ 


0.08+ 


0.06 


HD160641 


35.5 


0.45 


2.7+ 


4.5+ 


11.8 


500 


-7.3 





1.0 


0.3 


0.13 


0.11 


BD-9°4395 


25.1 


0.33 


2.5t 


4.4+ 


13.4 


400 


-7.9 


0.04+ 


1.3+ 


0.09+ 


0.14+ 


0.015+ 


BD+10°2179 


18.5 


0.06 


2.6t 


3.6+ 


12.3 


400 


-8.9 


0+ 


3.0+ 


0.08+ 


0.0005 


0.006 


HD144941 


27.0 


0.31 


3.9+ 


2.7+ 


8.4 


500 


-9.8 


1.4+ 


0.005+ 


0.003+ 


0.001 


0.01 + 


Sun* 


















0.21 


0.06 


0.066 


0.11 



+ taken from (or at least oriented at) previous photospheric analyses (cf. Table|5J 
t Solar metal abundances from Asplund et al. (2005) are given for comparison 



(i>d = 50 km s ) and for the photospheric spectrum in the optical 
(25kms _1 ). 

The UV spectrum of HD160641 shows wind-line profiles of 
the Si IV and C IV resonance doublets (see Fig.(6](, both embed- 
ded in a pronounced forest of iron lines. The P-Cygni profiles 
are not very wide, and fitted with a moderate wind velocity of 
Voc — 500 km s" 1 . A two-/3-law, where 40% of the velocity is 
attributed to a slow acceleration with /9a = 4, has improved the 
nice fit of the profile shape. 

By adjusting the mass-loss rate, the model matches both lines 
simultaneously. Again, the fit is not very sensitive to higher values 



because of saturation. Also for this star the "work ratio" is close 
to unity (1.15) and supports the plausibility of a radiation-driven 
wind. 

The observation also shows the Nv resonance doublet as a P- 
Cygni profile, while the model only predicts photospheric absorp- 
tions. Due to the relatively low stellar temperature, the ionization 
stage of N V is not populated in the stellar wind. Hence the observa- 
tion of that wind line is a hint that "superionization", known from 
the winds of massive O stars, is also effective in the winds studied 
here. 





1520 
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Figure 4. Normalized line spectram for BD+37° 1977 as in Fig.[3] zoomed 
at the UV resonance doublets of N V (top) and C IV (bottom). Stellar lines 
are marked in black, interstellar lines (i.s.) and IUE Resau marks (R) are 
marked in green. 



4.4 BD-9°4395 (V2209 Oph) 

BD— 9°4395 is a somewhat weird star, as it shows variations in to- 
tal light, radial velocity and line spectrum. These peculiarities have 
been attributed to non-radial oscillations and dynamical activities 
in a circums tellar shell. A detailed analysis of this star has been 
published bv lJefferv & Heberl dl992h . 

We adopt a slightly higher effective temperature (25.1 kK) 
than this previous study (22.7 kK) because otherwise our model 
predicts too strong lines from low ions like C II. Abundances have 
been determined in the previous photospheric analysis and can be 
adopted from there, leading to a good fit of the photospheric line 
spectrum (apart from the mentioned weird features) including the 
UV iron line forest (see Online Material, Fig.[3j- The SEP fits per - 
fectly when adopting the reddening law from Cardell i et al. I Gall) 
with Rv = 2.8 (which gives a marginally better match to the IUE 
data than the standard Rv = 3.1), and reveals a higher interstellar 
reddening than found in previous work. 

The stellar wind manifests itself only marginally in the spec- 
trum of BD— 9°4395. The C IV resonance doublet forms a narrow 
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Figure 5. Normalized line spectrum for BD+37°442 with the UV reso- 
nance doublets of N V (top) and C IV (bottom). Observations are in blue 
and green, while the best-fit model is red-dotted. 



P-Cygni profile, which can be fitted with a terminal wind velocity 
of 400 km s _1 and suitably adjusted mass-loss rate. However, even 
with a "two-beta-law" for the velocity field (40% with P2 = 4), 
which worked fine for HD 160641, the round shape of the absorp- 
tion feature is not reproduced, and the calculated P-Cygni emis- 
sion is weaker than observed. Because of the small wind velocity, 
we calculated the formal solution with a Doppler broadening of 
hd = 30 km s _1 in the UV and 15 km s _1 in the optical. 

The Si IV resonance doublet shows blue-shifted, asymmetric 
absorptions which are also not perfectly reproduced by the model. 
We attribute the remaining mismatch to the neglect of pressure 
broadening for resonance lines in the PoWR code. 

The final model for BD— 9°4395 has a "work ratio" of 1.06 
and is thus nearly hydrodynamically consistent. 

4.5 BD+10°2179 

This star is the coolest of our sample. Slightly discrepant val- 
ues have been published for its effective temperature (19.5 and 
16.9 kK, respectively, cf. Table|2](. Therefore we inspect the flux 
distribution and the photospheric line spectrum. The SED fit for 
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Figure 6. Normalized line spectrum for HD 160641 with the UV resonance 
doublets of Si IV (top) and C IV (bottom). Observations are in blue, while 
the best-fit model is red-dotted. 
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Figure 7. Normalized line spectrum for BD— 9°4395 with the UV reso- 
nance doublets of Si IV (top) and Civ (bottom). Observations are in blue, 
while the best-fit model is red-dotted. 



models with T» < 18 kK implies Eb-v = 0.00, which is un- 
likely and in conflict with the interstellar Lya absorption. However, 
T» > 18 kK makes C II line at 1335 A weaker than observed. Only 
with T„ = 19 kK, the Stark-broadened Hel line wings fit nicely 
(keeping logg = 2.6 from previous work). As a compromise we 
finally adopt T efi = 18.5 kK. 

The element abundances also differ between the two previous 
photos pheric analyses. Starting with the values from lPandev et al.l 
(2006), we found that silicon must be even less abundant to avoid 
too strong photospheric absorption features of the Si IV resonance 
doublet. Fo r carbon we prefer the slightly higher abundance from 
lHebeiUl983l) . 

Only weak indications of a stellar wind are found in the UV 
spectrum of BD+10°2179. The Si IV resonance doublet shows 
blue-shifted, asymmetric absorptions. When adopting a two-beta 
velocity field with Voo — 400 kms - and a suitable mass-loss rate, 
the model can roughly reproduce these profiles (Fig.[8}. But since 
the wind profile is not fully developed, v x is not accurately con- 
strained. Because of the weak wind, we calculated all non-iron lines 
in the formal solution with a small microturbulence of 15 km s _1 . 

From the C IV resonance doublet, only the red component 



seems to be reproduced by the same model. However, the feature 
apparently fitting the red C IV component is in fact due to iron. This 
is demonstrated in Fig.[8]by the thin, green dashed profile which is 
calculated with the iron forest switched off. Thus we must state 
that the model does not reproduce the C IV resonance doublet as 
a wind line at all. We found no escape from this discrepancy. A 
(much) larger mass-loss rate makes not only the Si IV resonance 
doublet too strong, but also transforms the C II line at 1335 A into 
a wind-line profile which is not observed. Reason for the missing 
C IV wind line is that models for such low T Q g predict that carbon 
recombines to C III and C II. We must conclude that also in this star 
some "superionization" keeps the ionization higher than predicted 
by our stationary models. As we will discuss below ((5), frictional 
heating in this very thin wind may provide a viable mechanism. 

Our "final model" and the mass-loss rate thus relies only 
on the Si IV resonance doublet, and should be taken with care. 
The "work ratio" is 0.41, i.e. not far from what is expected for a 
radiation-driven wind. 
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Figure 8. Normalized line spectrum for BD+10°2179 with the UV reso- 
nance doublets of Si IV and C IV. Observation is in blue, while the best-fit 
model is red-dotted. The thin-dashed green profiles are calculated without 
iron opacities in the formal integral. 



Figure 9. Normalized line spectrum for HD144941 with the UV resonance 
doublets of Si IV (top) and C IV (bottom). Observations are in blue, while 
the best-fit model is red-dotted. 



4.6 HD144941 

Previous photospheric analyses give quite discrepant effective tem- 
peratures for this star (see Table d . From a detail ed analysis of the 
optical line spectrum . lHarrison & Jeffervl d 19971) derived T c s = 
23.2 kK. In contrast, Jeff erv et alj feOOll) obtained 27.8 kK from 
fitting the UV flux distribution. We adopt T eS = 27.0 kK as 
a compromise, but take t he gravity and derived luminosity from 
Harrison & Jeffervl dl997b . Some of the pressure-broadened He I 
lines appear too broad when compared to our model (see Fig.|5]in 
the Online Material). 

The chemical abundances were als o determined from the pho - 
tospheric spectrum in previous work dHarrison & Jeffervl 1 1997b . 
Comparison with our model (see Fig.[5] in the Online Material) 
gives the impression that the traces of hydrogen might have been 
under-estimated. The formal integral is calculated for the non-iron 
lines with a microturbulence of vd — 30kms _1 in the UV and 
lSkms" 1 in the optical. 

The spectral signatures of the stellar wind are seen in the UV 
resonance doublets of Si IV and Civ (see Fig.[9j. In both cases, 



there are blueward extended absorption features accompanied by 
weak, red-shifted P-Cygni emissions. The reproduction of the wind 
profiles requires a terminal wind speed of about 500 kms - . The 
mass-loss rate of the final model (see Table|3j is a compromise, 
leaving the wind feature of the C IV doublet too weak. 

For both resonance doublets, the model shows too strong 
photospheric absorption features at the rest wavelengths shining 
through the wind. In the case of Si IV this discrepancy was so strong 
that we felt the need to reduce the silicon abundance against the 
value from the literature (from 25 to 10 ppm mass fraction), but the 
unshifted features are still stronger than observed. 

We also tested the low er effective temperature (23.2 kK) from 
lHarrison & Jeffervl i 19971) for the wind-line fit. With these mod- 
els, the problem with the too strong photospheric components al- 
most disappears. However, T e ff = 23.2 kK is too low to produce 
the observed the C IV wind line in our radiative-equilibrium mod- 
els. In this case, one would need to argue that C IV is due to some 
"superionization", as in the case of N V in HD 160641 and C IV in 
BD+10°2179 (see above). 

These ambiguities make the derived mass-loss rate especially 
uncertain. The work ratio for the final model is only 0.2, possibly 
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Figure 11. Mass-loss rate versus luminosity for the studied extreme he- 
lium stars (labelled red symbols), compared to various other groups of hot 
stars from high to low masses. Grey groups of symbols denote massive 
Wolf-Rayet stars (squares: WNE, diamonds: WNL). Also in grey are central 
stars of planetary nebulae of hydrogen-rich composition (CSPN [O], dots) 
or of Wolf-Rayet type (squares: [WCE], diamonds: [WCL]). The long blue 
dashed line gives the slope of M a L 1,5 , tight to the protopype O-type 
supergiant £ Puppis. 



indicating that the mass-loss rate is over- or the luminosity under- 
estimated. 



4.7 Wind variability 

Many extreme helium stars are known to be intrinsically variable 
in bot h light and radial velocity on timescales of one to several 
days dLvnas-Grav et aljl987l ; |jefferv & Heberl 19921) . Repeat high- 
resolution spectra of HD160641 and BP— 9°43 95 were obtained 
with IUE on timescales of days and years. I Jeff erv & Heberj 1 1992) 
found no evidence for variability of the wind lines in BD— 9°4395. 
FigurefTOlshows the mean profile of two resonance lines with strong 
wind profiles in boths stars. It also shows the standard deviation of 
the individual observations around the mean. While there is a hint 
of increased scatter in the blue wing of some lines, it is no greater 
than seen in other nearby regions of spectrum. It is concluded that 
there is little firm evidence for any variation in the wind line pro- 
files. In contrast, p hotospheric lines are known to show show strong 
profile variability jjefferv & Heberlll992l ; lwright et alj200r3) . 



5 CONCLUSION 

The spectra of six extreme helium stars under investigation are 
found to indicate weak mass loss. This evidence comes from reso- 
nance lines in the UV, which show P Cygni profiles or at least asym- 
metric absoption features. The hotter two of our programme stars, 
with T e g > 45 kK, exhibit wind profiles of the N IV and C IV res- 
onance doublet, while the cooler stars (T e g < 35 kK) show these 
effects in Si IV and C IV. 



By fitting these wind lines to model spectra, we have deter- 
mined the mass-loss rates. The result depends inversely on the 
abundance of the respective element, which we have mostly taken 
from previous analyses of the photospheric spectra. Thus one must 
realize that, apart from the uncertainties of the fit, a possible error 
in the element abundance directly translates into a wrong mass-loss 
rate. 

One should also be aware that the empirical mass-loss rate 
scales with the adopted distance to the star. The reason is that the 
resonance-line fits basically measure an optical depth which, in a 
homologuous scaling, depends on the density times the linear di- 
mension. The density depends on M /Ft 2 , and therefore the optical 
depth scales as M/R,. For models with the same optical depth, we 
thus obtain M oc R and, since L oc R 2 , a scaling M oc L 1 ^ 2 . 

Moreover, the fit relies on the ionization structure of the mod- 
els calculated in radiative equilibrium. When the ion under consid- 
eration is abundant, this should be reliable. However, in some of 
the cooler stars carbon is mainly in the stages C II and C III. In such 
cases the traces of the observed C IV ion can be largely enhanced 
by additional ionization processes, such as that caused by hydrody- 
namic shocks in the stellar wind which emit radiation in the hard- 
UV and X-ray range. Such "superionization" is well known from 
the winds of massive hot stars. Mass-loss rates cannot be reliably 
derived from "superionized" lines. 

Thus we warn that the mass-loss rates should be taken with 
care, although the error margin is impossible to quantify for the 
individual stars. 

As a general result, stellar winds appear as a persistent feature 
of hot stars, irrespective of their wide range of luminosities, masses 
and chemical composition. In Fig.QTJwe plot the empirical mass- 
loss rates M over the stellar luminosity L, both on a logarithmic 
scale. The theory of ra diation-driven winds in its most elementary 
form JCastoret all 19751) predicts a correlation of about M oc L , 
which roughly remains valid in later refinements of this prediction 
such as t he "modified wi nd momentum - luminosity relationship" 
(see e.g. IPuls et alj|2008h . Fig.QT] reveals that this relation, repre- 
sented by the blue dashed line normalized to the values for the mas- 
sive O-type supergiant £ Puppis, indeed marks a kind of lower limit 
for the mass loss from hot stars of all kinds. 

Our extreme helium stars lie close to this relation. Among the 
massive stars, the Wolf-Rayet types show considerably stronger 
mass loss. Similarly high M are found for the early subtypes 
(WNE) and the late subtypes (WNL), but the latter ones have on 
the average higher luminosities. 

Mass loss is also found from central stars of planetary nebu- 
lae (CSPN), which are hot stars of lower mass and less luminos- 
ity than the massive early-type stars. In their majority the CSPN 
are hydrogen rich (CSPN-[0]). Their mass-loss rates lie about one 
order of magnitude above the dashed M -L-relation in Fig.QT] Ex- 
tremely strong is the mass loss from the hydrogen deficient CSPN 
with Wolf-Rayet type spectra. According to detailed hydrodynamic 
modelling, the extremely strong mass loss from Wolf-Rayet stars is 
caused by their proxim ity to the Eddington limit , and driven by 
multiple line scattering (Grafener & Hamann 200^). 

The proximity to the Eddington limit can be measured by the 
ratio between the inward force by gravitation and the outward force 
by radiation pressure, V = g ra< j / <? gra v • Only accounting for the ra- 
diation pressure on free electrons by Thomson scattering, the Ed- 
dington r follows from the effective temperature and the surface 
gravity as 



1n -l5.l2 rp4 / 

10 q T ef f/g g 



(1) 
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Figure 10. The Si IV and C IV ultraviolet resonance doublets for two extreme helium stars known to be intrinsic variables as observed with IUE. The black 
histogram shows the mean obtained by combining all IUE HIRES and HST STIS images at these wavelengths. The red polyline shows the standard deviation 
about this mean. 



Table 4. Ratio between the inward force by gravitation and the outward 
force by radiation pressure, T = g r ad / Sgrav ■ 



Star 


r 


HD 160641 


0.55 


BD-9°4395 


0.32 


BD+37°1977 


0.17 


BD+37°442 


0.17 


BD+10°2179 


0.10 


HD144941 


0.03 



where T e g is in Kelvin and g in cgs-units. For our stars we may set 
the number of free electrons per atomic mass unit to q — 0.5 for 
fully ionized helium. 

To evaluate F, we must realize that the gravities g given 
in Table[3] are spectroscopically determined from the pressure- 
broadening of lines, and thus are effective gravities g = g g rav(l — 
r). Taking this into account, we obtain the V values compiled in 
Table|4] The correlation with the mass-loss rates is striking, espe- 
cially within the group of four stars with nearly the same lumi- 
nosity. The outstandingly high mass loss from HD 160641 is thus a 
consequence of its low gravity. 

This agreement with the expectations for radiation driven 
winds, also reflected by the plausible "work ratios" for the indi- 
vidual stars as mentioned in Sect.[4] may be considered an inde- 
pendent argument that the luminosities and distances are correct at 
least within an order of magnitude. Adopting a different luminos- 
ity would shift the empirical mass-loss rate parallel to the thin line 
with slope 1/2 indicated in the lower-right comer of Fig.ll II 

It has been theoretically predicted that in very thin winds 
the radiation pressure will accelerate only the metal ions, while 
the bulk matter of helium stays inert. This ion decoupling will 
lead to frictional heating, or perhaps eve n to instabilities that com- 
pletely disrupt the smooth ste llar wind iSpringmann & Pauldrachl 
1 19921 ; iKrticka & Kubatll20o"lh . Among our programme stars, the 



thinnest wind is encountered at BD+10°2179. According to our 
final model, the number density at one stellar radius from the pho- 
tosphere is about 10 8 atoms per cubic centimeter. The mass- 
loss rate and radius of that star are almost identical to those of 
r Sco, for which ion decoupling has been predicted. Hence we may 
speculate that the "superionization" invoked for C IV (Sect. 14.5^ is 
caused by frictional heating from ion decoupling. The same may 
also hold for HD 144941, the star with the second-thinnest wind 
(log n a tom/cm~ 3 = 8.4 at r = 2 R t ). The consistent wind-line 
fits for the other programme stars provide evidence that ion decou- 
pling does not take place in their winds. 
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Figure 1. BD+37°442. Top panel: Spectral energy distribution (SED) for from the UV to the near IR. The observations comprise low-resolution IUE data 
(thin blue line), visual photometry (blue blocks with magnitudes) and near-IR photometry from 2MASS. The model continuum (see Table[3]for the parameters) 
is shown by the straight red line, in parts of the spectrum augmented by the synthetic line spectrum. Middle panel: Line spectrum in the UV; the observation 
(blue line) has been divided by the model continuum for normalization. The model spectrum is shown by the red-dotted line. Bottom panel: Same for the 
optical range. This observation has been normalized to the continuum "by eye". 
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